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' The topological susceptibility of QCD vacuum is studied in the framework of a 

covariant chiral quark model with non-local quark-quark interaction. The relation 
of the first moment of topological susceptibility x'(0) an d the 'spin crisis' problem 
is briefly discussed. It is shown, in particular, that one always gets the inequality 
X'(0) > X'ozi- 
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It is well known that due to U a (1) axial Adler-Bell-Jackiw anomaly the iso-singlet 
t^. . axial-vector current 

J$s =H 9/7/i759/ (1) 
is not conserved even in the chiral limit, and its divergence equals 

^ vS? (*) = 2N fQ* (*) . ( 2 ) 

Ph. where 

Q 5 (x) = (aj8n)G«„(x)G« v (x) (3) 
is the topological charge density. The correlator of singlet currents is defined as 

X ' 

= */ d 4 xe^(0\T{j$(x)J%(0y}\0) = 
= (q,qu-9,uq 2 )n%(Q 2 ) + q,q u U^ L (Q 2 ). (4) 

In the chiral limit the longitudinal part of the correlator defines the topological suscepti- 
bility, i.e. the correlator of the topological charge densities, Q$ (x), 

x(Q 2 )=i J d 4 xe^(0|T{Q 5 (x)Q 5 (0)}|0>, (5) 

with the relation (see, e.g., [1]) 



n£'°(Q 2 ) = ^4(Q 2 ). (6) 



At high Q the operator product expansion (OPE) predicts [2] 

a„ /a. /„„ \2\ 



4« 2 --) = -Ts(?( G -)) +0 ' Q " 2 )' < 7 > 
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where the perturbative contribution has been subtracted. 

At low Q 2 , x (Q 2 ) is represented as a sum of contributions coming purely from QCD 
and from (tt, rj) — mesonic resonances [1] 
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m u m d 
m u + m d 
f 2 



(uu) - X '(0)Q 2 



(8) 



m u - m d 
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4 [\m u + m d J Q 2 + ml 3Q 2 + m 2 _ 



+ 0(Q 4 ), 



where the first term has been found in [3] and its chiral limit follows the Crewther theorem 
[4] maintaining that % (0) = in any theory where at least one massless quark exists. 
The estimates of x'(0) existing in the literature are rather controversial: 



X'(0) = (48 ± 6 MeV) 2 [5], x'(0) = (26 ± 4 MeV) 5 



(9) 



Both estimates were found within the QCD sum rules method. These values of the first 
moment of topological susceptibility have to be compared with the value obtained in the 
Okubo-Zweig-Iizuka (OZI) case, the case free of axial anomaly which is 



Xozi (0) 



J TV 

2N f 



(39 MeV) 2 



The principal point is that smallness of x'(0) is the base for the one of the mechanisms 
explaining "proton spin crisis" problem [6]. Indeed, within this approach it is assumed 
that the flavor singlet axial charge ao {Q 2 ) is proportional to the product of the first 
moment of the QCD topological susceptibility taken at scale Q 2 and an RG-invariant 
coupling of 'OZI Goldstone boson' with nuclon 



a (Q 2 ) 



2m 



N 



(10) 



This mechanism has been, however, criticized in [8]. All this makes important further 
model estimates of x'(0). 

Within the chiral quark model 2 [9] based on the non-local structure of instanton QCD 
vacuum [11] the full iso-singlet axial- vector vertex becomes [10] 



T%(k,q,k' = k + q) 



7^ - (k + k ) lt ±- 



r.2 



2^M(k')M(k)-- 



k' 2 - k 2 
G' 1 - GJ PP (q 2 ) 



(11) 



G'J PP (q 2 



75- 



where M(k) is dynamical, momentum dependent quark mass, G and G' are 4-quark 
couplings in iso-triplet and iso-singlet channels, correspondingly, and 



Jp P {q 2 )S ab 



d 4 k 



M 2 J (27r)' 



M (k) M(k + q) Tr S(k)^T a S (k + q) -y 5 T b 



;i2) 



2 The explicit calculations below are performed in SU (2) sector of the model. 
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In (12) the (inverse) quark propagator is S 1 {p) = ft — M{p). Because of axial anomaly 
the singlet current does not contain massless pole, since as q 2 — > one has: 



1 - GJ PP (q 2 ) _ ft 
-q 2 M\ ' 



(13) 



where f n is pion weak decay constant and M q = M(0). The cancellation of the massless 
pole occurs with help of the gap equation. Instead, the current develops a pole at the r( — 
meson mass 3 , 1 — G'J PP (q 2 = —m 2 ,) = 0, thus solving the C/a(1) problem. The vertex 
(11) satisfies the anomalous Ward-Takahashi identity: 



m , , 2JM(k')M(k) ( 

q,Y { % (k, q,k' = k + q) = 75^ (k') + S? 1 (k) 7s + 75 Q IJpp ^ [ X G 



(14) 



where the last term is due to the anomaly. Thus, the QCD pseudoscalar gluonium operator 
is interpolated by the pseudoscalar effective quark field operator with coefficient expressed 
in terms of dynamical quark mass. This is a consequence of the fact that in the effective 
quark model the connection between quark and integrated gluon degrees of freedom is 
fixed by the gap equation. 

For completeness we display the vertex corresponding to the conserved iso-triplet axial- 
vector current 



r« 5 (k,q,k' = k + q)=T a 



7m - % 



M(k') + M(k) 



(15) 



k' 2 - k 2 

satisfying the axial Ward-Takahashi identity 

q.T^k, q, k') = l5 S^ (k 1 ) T a + T^ 1 (k) 7s . 



75 



(16) 



The axial- vector vertex (15) has a kinematical pole at q 2 = 0, a property that follows 
from the spontaneous breaking of the chiral symmetry in the limit of massless u and d 
quarks. Evidently, this pole corresponds to the massless Goldstone pion. 

The quark matrix elements of currents corresponding to vertices (11) and (15) can be 
expressed in terms of real form factors 



p's' 



4°5 3) (0)| ps) = Us' (PO [ 7 ,7 5 GS°' 3) (q 2 ) - q^G^ (q 2 )} u s (p) , (17) 



where = (l,r 3 /2), u s (p) are spinor solutions of the Dirac equation for free quarks, 

and the currents are defined as 



4°5 3) (l) = I 7^4^ ( k ) q,k' = k + q)i, (k + q), 



with ip (k) being the solutions of the Dirac equation 

k-M(k)]ip(k) =0. 



(18) 



(19) 



3 See previous footnote. 
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By using the Dirac equation one gets 

8"42 = 0, 9M<S (g 2 ) = i ig^ (g2) / -j^2y/M(k>)M(kW 75 ^ {k + g ) 

Comparison with (17) leads to the relations for form factors (taken in the local limit 

M(k) « M,) 

GS 3) (g 2 ) = 1, (g 2 ) = 2M 9 /g 2 , gS 0) (0) = 1, Gf (0) = 

resembling the results for a model of free massive quarks. 

Full model calculations lead to the following expression for the topological suscepti- 
bility [10] 



-(2iV» 2 X (Q 2 ) 



Q 2 J, A (Q 2 ) 



G' 



+ M 2 J pp (Q 2 )[1-—J ap (Q 



_ G'J AP (Q 2 ) 1 

Ml + 1 - G'J PP (Q 2 ) 
~GJ AP {Q 2 ) G-G' 



(20) 



M 2 

q 



G[1-G'J PP (Q 2 )] 



+ 



G 
H 2 

q L 



AN c N f 



d A k M(k) 
(2vr) 4 D (k) 



M (k) - \JM (k + Q) M (k) 



where D(k) = k 2 + M 2 (k) and the integrals J AP (q 2 ) and J^^q 2 ) are defined by 
J AP (q 2 ) = W c N f J (I + q) M (I), 



J*A (q 2 ) S ab = % I 4r^ Tr \ S ( k K a ( k ' 9, k + q)S (k + q) Yl (k + q, k) 
v ' q J (27r) l 

At large Q 2 one obtains the power-like behavior consistent with the OPE prediction (7), 
namely 



(2N f ) 2 X (Q 2 ^oc) 



2N f M 2 ( 1 _<?\ 



G 



G 



At zero momentum the topological susceptibility is zero 

x(o) = o, 



(21) 



(22) 



in accordance with the Crewther theorem. For the first moment of the topological sus- 
ceptibility we obtain [10] 



(23) 



If the OZI rule were exact in the flavor singlet channel and there were no anomaly, then 
one would have G' — G and x' (0) = Xozi (0)- But actually one has strong attraction in 
the iso-triplet channel and strong repulsion due to the anomaly in the iso-singlet channel 
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that means that one has always inequality G' < G. The second, negative term in (23) 
is numerically suppressed with respect to the first, positive term, J^p(O)/ f% f 



-0.24. 



Thus, from the existence of the anomaly we always have inequality 



x' (0) > X'ozi (0) , (24) 

and it is impossible to get anomalously small x' (0)- At this point we also mention other 
alternative approaches to the spin crisis problem based on screening of topological charge 
in the QCD vacuum [12, 13] (see for review, e.g. [14]). 

The constants G and G' are fixed with the help of the meson spectrum. Approximately 
one has G' ~ 0.1 G. As a profile for the dynamical quark mass we take a Gaussian form 

(25) 



M{u) = M 9 exp (-2m/ A : 



with the model parameters M q = 0.3 GeV, A = 1.085 GeV. Then the estimate for the 
first moment of the topological susceptibility is [10] 

X '(0) = (50 MeV) 2 . (26) 

To get the above result we have taken Nf = 3 in Eq. (23). We can see that the model 
gives the value of x'(0) which is close to the estimate of Ref. [5]. The influence of the 
current quark masses on x'(0) is expected to be small and the contribution of it- and 
//-mesons may be found from Eq. (8) 



XtT' 



(0) » (28 MeV) 5 



x (d 1 )(ld'GeV t ) 



and thus Xtoti®) ~ (57 MeV) 2 for the total result. 

The model prediction for the topological susceptibility is shown in Fig. 1. In the region 
of small and intermediate momenta our result is quantitatively close to the prediction of 
the QCD sum rules with the instanton effects included [1]. 

In the present talk we analyzed the cor- 
relation function of the singlet axial-vector 
currents within an effective non-local chiral 
quark model. By considering this correla- 
tor the topological susceptibility was found 
as a function of the Euclidean-momentum 
and its first moment was estimated. We 
demonstrated that in realistic situation one 
always gets the inequality x'(0) > xbzi(®)i 
thus discarding the mechanism explaining 
the 'spin crisis' based on anomalous small- 
ness of x'(0)- m addition, the fulfillment of 
the Crewther theorem was demonstrated. 
It would be interesting to verify the pre- 
dictions given in Fig. 1 by modern lattice 
simulations. 
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partially supported by grants RFBR (nos. 02-02-16194, 03-02-17291) and INTAS-00-00- 
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Figure 1. Topological susceptibility versus Q 2 
predicted by the model with G' = 0.1 G, Eq. 
(20). 



5 



References 

[1] B.L. Ioffe, arXiv: hep-ph/9811217; B.L. Ioffe, A.V. Samsonov, Phys. Atom. Nucl. 63 
1448 (2000) [Yad. Fiz. 63 1527 (2000)]. 

[2] V.A. Novikov, MA. Shifman, A.I. Vainshtein and V.I. Zakharov, Phys. Lett. B 86 
347 (1979). 

[3] G. Veneziano, Nucl. Phys. B 15 2139 (1979); P. Di Vecchia, G. Veneziano, Nucl. 
Phys. B 171 253 (1980) . 

[4] R.J. Crewther,Phys. Lett. B 70 349 (1977). 

[5] B.L. Ioffe, A.G. Oganesian, Phys. Rev. D 57 6590 (1998). 

[6] S. Narison, G.M. Shore and G. Veneziano, Nucl. Phys. B 546 235 (1999). 

[7] G.M. Shore, G. Veneziano, Phys. Lett. B 244 75 (1990). 

[8] M. Anselmino, A. Efremov and E. Leader, Phys. Rept. 261, 1 (1995). 

[9] A. E. Dorokhov and L. Tomio, Phys. Rev. D 62 014016 (2000); I. V. Anikin, 
A. E. Dorokhov and L. Tomio, Phys. Part. Nucl. 31 509 (2000) [Fiz. Elem. Chast. 
Atom. Yadra 31 1023 (2000)]. 

[10] A. E. Dorokhov and W. Broniowski, Eur. Phys. J. C 32 (2003) 79. 

[11] A.E. Dorokhov, S.V. Esaibegian, S.V. Mikhailov, Phys. Rev. D 56 4062 (1997); A.E. 
Dorokhov, S.V. Esaibegian, A.E. Maximov, S.V. Mikhailov, Eur. Phys. J. C 13 331 
(2000) 

[12] S. Forte, Phys. Lett. B 224, 189 (1989). 

[13] A.E. Dorokhov, N.I. Kochelev, Mod. Phys. Lett. A 5 55 (1990). 

[14] A.E. Dorokhov, N.I. Kochelev, YA. Zubov, Int. J. Mod. Phys. A 8 (1993) 603. 



6 



